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Abstract. We develop the analysis of finite element approximations of im- 
plicit power-law-like models for viscous incompressible fluids. The Cauchy 
stress and the symmetric part of the velocity gradient in the class of models 
under consideration are related by a, possibly multi-valued, maximal monotone 
r-graph, with 1 < r < oo. Using a variety of weak compactness techniques, 
including Chacon's biting lemma and Young measures, we show that a subse- 
quence of the sequence of finite element solutions converges to a weak solution 
of the problem as the finite element discretization parameter h tends to 0. A 
key new technical tool in our analysis is a finite element counterpart of the 
Acerbi-Fusco Lipschitz truncation of Sobolev functions. 



1. Introduction 

Most physical models describing fluid flow rely on the assumption that the 
Cauchy stress is an explicit function of the symmetric part of the velocity gra- 
dient of the fluid. This assumption leads to the Navier-Stokes equations and its 
nonlinear generalizations, such as various electrorheological flow models; see, e.g., 
|Lad69l ILio69l IRiizOO| . It is known however that the framework of classical con- 
tinuum mechanics, built upon the notions of current and reference configuration 
and an explicit constitutive equation for the Cauchy stress, is too narrow to en- 
able one to model inelastic behavior of solid-like materials or viscoelastic properties 
of materials. Our starting point in this paper is therefore a generalization of the 
classical framework of continuum mechanics, called the implicit constitutive theory, 
which was proposed recently in a series of papers by Rajagopal; see, for example, 
|Raj03[ |Raj06| . The underlying principle of the implicit constitutive theory in the 
context of viscous flows is the following: instead of demanding that the Cauchy 
stress is an explicit function of the symmetric part of the velocity gradient, one may 
allow an implicit and not necessarily continuous relationship between these quanti- 
ties. The resulting general theory therefore admits fluid flow models with implicit 
and possibly discontinuous power-law- like rheology; see, |Mal08[ [MalO?! . Very re- 
cently a rigorous mathematical existence theory was devel oped for the se models by 
Bulfcek, Gwiaz da, Malek, a nd Swierczewska-Gwiazda in |BGMSG09| . Motivated 
by the ideas in BGMSGOO] . we consider the construction of flnite element approx- 



imations of implicit constitutive models for incompressible fluids and we develop 
the convergence theory of these numerical methods by exploiting a range of weak 
compactness arguments. 
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Let f7 c M'', fi e N, be a bounded open Lipschitz domain with polyhedral bound- 
ary. For re (1, oo) we define r' e (1, oo) by the equality ^ + -p = 1 and we set 

f := min{r',r72}, where r* := i ^ if r < d, 

1 00 otherwise. 

We refer the reader to the first paragraph in Section 12.11 for the definitions of the 
various function spaces used in the paper and for a list of our notational conventions. 

Problem. For / e W-^'''\ny find {u,p, S) e Wa'^^ifl)"^ x L~^{n) x W' {nY""^ such 
that 

<X\Y{u®u + pl- S) = f in V'{n), 

divM = in V'{VL), (1.2) 
{Du{x), S{x)) 6 A{x) for almost every x e ft. 

The symmetric part of the gradient is defined by Du := ^(Vm + (Vm)"^). As is 
implied by the choice of the solution space WQ'^{n)'^ for the velocity, a homogenous 
Dirichlet boundary condition is assumed here for u. 

The implicit law, which relates the shear rate Du to the shear stress S, is given 
by an nonhomogeneous maximal monotone r-graph A : x ^-^ A{x). In particular, 
we assume that the following properties hold for almost every x e ft: 
(Al) (0,0)e^(a;); 
(A2) For aU {Si, en), {S2, (T2) e A{x), 

{{Ti — (T2) : {Si — S2) 5= {-^{x) is a monotone graph), 
and if 7^ ^2 and cri ¥= <T2, then the inequality is strict; 
(A3) If (<5, cr) e Rfy^„f x MfyX„f and 

{a - cr) : {S - S) ;^ for aU (5,0-) e ^(x), 
then {S,(t) e A{x) (i.e., A{x) is a maximal monotone graph); 

(A4) There exists a constant c > and a nonnegative m e L^{n), such that for all 
{S, cr) e A{x) we have 

cr : S ^ c{\S\^ + \cr\^ ) — m{x) (i.e., A{x) is an r-graph); 

(A5) (i) For aU S e Rf^^f := {C e R'*'^'' ^ = C^}, the set 

{<TeMf;„^(J,^)e^(a;)} 

is closed; 

(ii) For any closed C cz R^ynu the set 

{{x,S) e Q, X R^y^ : there exists cr e C, such that {S, cr) e A{x)} 

is measurable relative to the smallest cr-algebra £(f2)(g)*8(Rfjf,j'f) of the 
product of the cr-algebra £(il) of Lebesgue measurable subsets of fl 
and all Borel subsets *B(R,^y^nf ) of W^^^. 

The class of fluids described by (|1.2p is very general and includes not only New- 
tonian (Navier-Stokes) fluids {S = 2^if,Dv with ^^f, being a positive constant), but 
also standard power-law fluid models, where S = 2^^\DvY~'^ Dv, 1 < r < 00, 
and their generalizations {S = 2jl{\Dv\'^) Dv), stress power-law fluid flow models 
and their generalizations of the form Dv = a{\S\'^) S, fiuids with the viscosity de- 
pending on the shear rate and the shear stress S = 2fl{\Dv\'^, \S\'^) Dv, as well as 
activated fiuids, such as Bingham and Herschel-Bulkley fluids. For further details 
concerning the physical background of the implicit constitutive theory we refer the 
reader to the papers by Rajagopal and Rajagopal & Srinivasa |Raj03[ [Raj06[ |RS08| . 
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and the introductory sec tions of Buhcek, Gwiazda, Malek & Swierczewska-Gwiazda 
|BGMSG09[ fBGMSGllI and Buhcek. Malek & Suh |BMS12j . 



It is proved in |BGMSG09 that under the assumption r > problem (|1.2p has 



a weak solution. The proof in |BGMSG09j uses a sequence of approximation spaces 
spanned by finite subsets of a Schauder basis of an infinite-dimensional subspace of a 
Sobolev space, consisting of exactly divergence- free functions. Since such a Schauder 
basis is not explicitly available for computational purposes, here, instead, we shall 
approximate p.2p from two classes of inf-sup stable pairs of finite element spaces. 
The first class contains velocity-pressure space-pairs that do not lead to exactly 
divergence-free velocity approximations. For finite element spaces of this kind our 
convergence result is restricted to r > -j^. In the case of exactly divergence- free 
finite element spaces for the velocity we show that the resulting (sub)sequence of 
finite element approximatio ns converges to a weak solution of the problem for any 
r > J^, as in |BGMSG09| . 

The paper is structured as follows. In Section [2] we introduce the necessary 
analytical tools, including Young measures and Chacon's biting lemma. In Section 
ini we define the finite element approximation of the problem with both discretely 
divergence-free and exactly divergence-free finite element spaces for the velocity. A 
key technical tool in our analysis is a new discrete Lipschitz truncation technique, 
which can be seen as the finite element counterpart of the Lipschitz truncation of 
Sobolev functions discovered by Acerbi and Fusco fAF88' and further refined by 
Diening, Malek, and Steinhauer |DMS08| : see also [DHHRll, BDF12[|BDST2] . The 
central result of the paper is stated in Section |4j in Theorem [181 ^^nd concerns the 
convergence of the finite element approximations constructed in Section [31 

2. Preliminaries 

In this section we recall some known results and mathematical tools from the 
literature. We shall first introduce basic notations and recall some well-known 
properties of Lebesgue and Sobolev spaces. We shall then discuss the approximation 
of an x-dependent r-graph by a sequence of regular single- valued tensor fields using 
a graph-mollification technique by Francfort, Murat and Tartar |FMT04| (which 
the authors of |FMT04| attribute to Dal Maso). We close the section by recalling a 
generalization from |Gwi05[rGZG07| of the so-called fundamental theorem on Young 
measures; cf. |Bal89) . 

2.1. Analytical framework. Let C{0.)'^ be the space of vector- valued continuous 
functions on Cl and let Co(IKfy^nf) denote the space of continuous functions with 
compact support in Kg^j^. For a measurable subset w cz M'*, we denote the classical 
spaces of Lebesgue and vector-valued Sobolev functions by (L*(a;),|| • \\s-uj) and 
{W^-'{ujy,\\ ■ \\i,s-u:), s e [l,oo], respectively. Let V{uj) := C^(a;)'' be the set 
of infinitely many times differentiable functions with compact support in lo; we 
denote by f '(w) the corresponding dual space, consisting of distributions on w. For 
s e [1, oo), denote by W^'^iuj^ the closure of in W^'''{uj)'^ and let W^;^i^{ujy := 
{v e Wq'^^uj)''' : div v = 0}. The case s = co has to be treated differently. We define 

and 

Moreover, we denote the space of functions in L'^{uj) with zero integral mean by 
For s,s' 6 (1, <») with ^ + p- = 1, L"^' {fl) and Lq (ft) are the dual spaces 
of L^n) and mn), respectively. The dual of '"(17)'' is denoted by W-^-^'iHy. 
For w = 57 we omit the domain in our notation for norms; e. g., we write \\-\\s instead 
of II • \\s,n- 
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Inf-sup condition. The inf-sup condition lias a crucial role in the analysis of in- 
compressible flow problems. It states that, for s,s' e (l,oo) with i + = 1, there 
exists a constant as > such that 

sup ^^^pl^ ^ as\\q\\s' for all q e (n); (2.1) 
compare with |DRS10[ |Bog79| . 

Korn's inequality. According to (|1.2p the maximal monotone graph defined in 



(Al) (A5) provides control only over the symmetric part of the velocity gradient. 



Korn's inequality states that this is already enough to control the norm of a Sobolev 
function; i.e., for s e (l,oo), there exists 7s > such that 

7s||-y||i,s \\Dv\\s for aU v e W^ '{n)''; (2.2) 
compare, for example, with |DRS10| . 

2.2. Approximation of maximal monotone r-graphs. In general an x-depen- 
dent maximal monotone r-graph A satisfying (Al) (A5) cannot be represented 
in an explicit fashion. However, based on a regularized measurable selection, it 
can be approximated by a regular single-valued monotone tensor field. Following 
|FMT04| . there exists a mapping S* : n x R'^^^^ R^^^ (a selection) such that, for 
all 5 e R^ym: (^^ S*{x, d)) e A{x) for almost every xen and 
(al) S* is measurable with respect to the product cr-algebra £(0) (g) 53(Msy^); 
(a2) For almost all x e the domain of S* is Rgymi 
(a3) S* is monotone, i. e., for every ^i, ^2 e ^tym almost all x e fl, 

{S*{x, S,) - S^x, S2)) : (^1 - 62) > 0; (2.3) 

(a4) For almost all x e il and all d e R^y^ the following growth and coercivity 
conditions hold: 

|S'*(a;,5)| s;ci|^r"^ + fc(a;) and S*{x,S) : d ^ C2\S\'^ - m{x), (2.4) 

where ci, C2 > 0, and k e (fl) and m e L^{fl) are nonnegative functions. 
Let 77 G Co(IRgy^j^) be a radially symmetric nonnegative function with support in 
the unit ball Bi{0) cz R^^^ and /gdxd 77 dC = 1. For n e N we then set 77" (C) = 
n'^ r]{n<^) and define 

5"(x, S) := (5* * 77")(x, S)= f S*(x, C) 7?"(^ - C) dC = / S*(^, C) dfi^- 

(2.5) 

Here, thanks to the equality Jjg^dxd -q" = 1 and the nonnegativity of 77. for each 

6 6 Mgy^Jn' d^s f7"(^ ~ C)dC defines a probability measure that is absolutely 

continuous with respect to the Lebesgue measure, with density 77" (5 — {■)). 

We reca ll the following properties of the matrix function S" from |GMS07[ 
lBGMSG09[r(;Z(;n7| . 

Lemma 1. The x-dependent matrix function S'", defined in (|2.5p . satisfies 

{S^{x, Si) - S'\x, S2)) : (5i - ^2) > for all 6^,52 G K^,"„t 

Moreover, there exist constants Ci,C2 > and nonnegative functions ifi e L^(Q,), 
k e (f2) such that, uniformly in n e N, we have 

\S'\x, S)\^h \Sr^ + kix) for all S e Rf;„t 

S"{x, 5):S^C2 \d\'' - rh{x) for all d e R^^^. 
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Remark 2. The selection S* enters in the definition of the finite element method 
in the form of through the Galerkin ansatz; compare with Section \3.4\ below. The 
natural question is then how one can gain access to such a selection. In fact, in 
most physical models it appears that the selection S* is given and A{x) is defined as 
the maximal monotone graph containing the set {{D, S*{D)) : D e M^^^}; compare 
with |BGMSG09[rGMS07| and the references therein. 



2.3. Weak convergence tools. The result in |Gwi05[ IGZG07| extends |Bal89) 
from limits of single distributed measures to limits of general probability mea- 
sures. To this end we need to introduce some standard notation from the theory 
of Young measures. We denote by OT(Rfy^^) the space of bounded Radon mea- 
sures. We call ji : VL ^ 931(Msjf„f), x i-^ fix, weak-* measurable if the mapping 
X 1-^ J^dx d h{C) dfixiC) is measurable for all h e Co{M.^y^^). The associated non- 
negative measure is defined by \fix \ (C) := fi^ + f^x ^ the Jordan decomposition 
IJLx = IJ't ~ t^Z i'^to t'^o bounded non- negative measures , ^~ . By means of the 
norm ||Mlli=o(n;OT(R?/„^)) '= esssup^^n /R^°i;i„^ d [/x^;! the space L'^{^;dn{W^^^)) of es- 
sentially bounded, weak-* measurable functions turns into a Banach space with 
separable predual L^(r2, Co(Mfjf^)). The support of a non-negative measure is de- 
fined to be the largest closed subset of Mfy^„f for which every open neighborhood of 
every point of the set has positive measure and supp /i^; := supp/ij u suppyU^. 

Theorem 3 (Young measures). Let Q. he an open and hounded suhset ofW^. Sup- 
pose that {ly^j^fii cz L^,{^l;dJl(S.'^y^^)) is such that v^. is a prohahility measure on 
l^sym fof 0,11 j G N and a.e. x e il. Assume that converges to v in the weak-* 
topology of U^ift; MiR^^^^D) for some v G L^(0; SH(Rf^>=„f)). 

Suppose further that the sequence {i'-'}jEN satisfies the tightness condition 

lim sup I {a; e 17 : supp iyi.\BR{0) 0} I ^ 0, 

where Bji{0) denotes the hall in Mf^„'f with center 6 Mf^„j and radius R > 0. 
Then, the following statements hold: 

(i) Vx is a prohahility measure, i. e., \\vx\\rfy^^^dxd^ = J^^dxd dlvxl = 1 o,.e. in fi; 
(it) for every h e i=^(r!; ^(Mf;^)), 

h{xX)di'i{C)—* [ h{xX)dvx{C) weak-* in L'^^ifl); 

d fmdxd 

(Hi) for every measurable subset uj cz ^} and for every Caratheodory function h 
such that 



lim sup / / \h{xX)\ dui{C)AC = Q (2.6) 

fl-*CO Jul J{C,^«.i^^: \h{xX)\>R} 

we have that 

h{x , C,) dvi.{C,) ^ / h{xX)'A^x{C) weakly in L^{uj). 

dx d I -mdx d 

sym ^ l^sym 

Lemma 4 (Chacon's biting lemma). Let il he a bounded domain in M'' and let 
{w"}„5=N be a hounded sequence in L^{fl). Then, there exists a nonincreasing se- 
quence of measurable subsets Ej c 17 with \Ej \ as j —> co, such that {v"}ni=N is 
precompact in the weak topology of L^(Q,\Ej) , for each j e N. 

3. Finite Element Approximation 

This section is concerned with approximating problem p.2p by a finite element 
method. To this end we introduce a general framework covering inf-sup stable 
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Stokes elements, which are discretely divergence-free, as well as exactly divergence- 
free finite elements for the velocity. These two classes of velocity elements require 
different treatment of the convection term. The discussion of these, including rep- 
resentative examples of velocity-pressure pairs of finite element spaces from each 
class, is the subject of ti3.2l and i j3.3l The finite element approximation of p.2p is 
stated in iJ3.4l We close with a new Lipschitz truncation method for finite element 
spaces, which plays a crucial role in the proof of our main result. Theorem [T51 

3.1. Finite element spaces. We consider a family {V",Q"}„eN c W^'"'\n)'^ 
xL"(J7) of pairs of conforming finite-dimensional subspaces of W^'^^fl)'^ x 
To be more precise, let {Gn}neN be a sequence of shape-regular partitions of H, i. e.. a 
regular finite element partition of Cl satisfying the following structural assumptions. 

• Ajfine equivalence: For every element E 6 C/„, n e N, there exists an invertible 
affine mapping 

Fe : E ^ E, 

where E is the closed standard reference d-simplex or the closed standard unit 
cube in R'^. 

• Shape-regularity: For any element E 6 C/„, n 6 N, the ratio of its diameter to the 
diameter of the largest inscribed ball is bounded, uniformly with respect to all 
partitions Qn, n e N. 

For a given partition Qn, n e N, and certain subspaces V c C{ft)'^ and Q c: 
the finite element spaces are then given by 

V" = V(a„) := {y e V : V^EoF^'e Py, E e g„ and V\,.n = o} , (3.1a) 

Q" = QiGn) := {q e Q : Q^e o F^' e F^j, E e e„} , (3.1b) 

where Py ^ W^''^''{EY and Pq c L'^{E) are finite-dimensional subspaces, with 
dimPy = ^ and dimPy = j, respectively, for some £,jeN. Note that Q" c L=^(fi) 
and since V" c C{n)'^ it follows that V" c: W^''^'{ny. Each of the above spaces is 
assumed to have a finite and locally supported basis; e. g. for the discrete pressure 
space this means that for n e N there exists Nn e N such that 

Q" = span{gi\...,g^J 

and for each basis function Q", i = 1, . . . , Nn, we have that if there exists E e Qn 
with # on E, then 

supp Q] cz\J{E' eQn\E' nE^ 0} =: il^. 

We introduce the subspace of V^;^ of discretely divergence-free functions by 

V^i, := {y 6 V" : <div V, Q}^ = for all Q e Q"} 

and we define 

Qo := {QeQ" : / Qdx = 0}. 
Jn 

Throughout the paper we assume that all pairs of velocity-pressure finite element 
spaces possess the following properties. 

Assumption 5 (Approximability). For all s e [l,co), 

inf \\v — V\\is^O for all V e Wn'^iil)^, as n ^ oo; and 

inf \\q — Q\\s for all q 6 L''*(fi), as n ^ oo. 

Assumption 6 (Projector njjj^). For each n e N there exists a linear projection 
operator n^;^ : Wo'\ny -> V" such that, for all s e (1, oo). 
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• li^iv preserves divergence in the dual of Q"; i. e., for any v e W^'^ijl)'^ we have 

<div v,Q)^= <div n^wv, Q)n for ah g e Q". 

• IIj;^ is locally VF^'^-stable; i. e., there exists ci > 0, independent of n, such that 

/ |VnSivW| dx ci / \Vv\ dx for all v e W^'i^)"^ and all E e g„. (3.2) 
Je Jq.^ 

Here we have used the notation j-^ ■ dx := ■ dx for the integral mean- value 

over a measurable set B ^ Mf^, \B\ ^ 0. 

It was shown in |BBDR10i IDR07| that the local W'l'i-stability of H^j^ implies 
its local and global VF^'^-stability, s e [l,oo]. In fact, by noting that the power 
function i i-^ is convex for s e [1, <»), we obtain by the equivalence of norms on 
finite-dimensional spaces and standard scaling arguments that 

\vn2,M^)\ *s ||vnji>|U,x,(s) *s c/ |vnSi,^| dx 

Je 

ikcf \Vv\dx iici j- \Vv\'' dx 

where we have used Jensen's inequality in the last step; recall that |r2^| ^ c \E\ with 
a constant depending solely on the shape-regularity of Qn- Raising this inequality 
to the s-tli power and integrating over E yields 

ivn^iv-yT dx^c [ \Vv\' dx. 

Summing over all elements E e Qn, accounting for the locally finite overlap of 
patches and noting Friedrichs' inequality, yields, for any s e [1, cx)), that 

||nji>||i,, s;c,||t;||i,, forall-(;6<-*(r!)^ (3.3) 

with a constant Cg > independent of n e N. Note that for s = cX) the inequality 
p.3p follows from an obvious modification of the argument above. 

Hence, by invoking the approximation properties of the sequence of finite element 
spaces for the velocity, stated in Assumption [5l we obtain that 

\\v - U'^i^vWi^, for aU v e W^'^inf, as n ^ cX) and s e [!,<»). (3.4) 

Assumption 7 (Projector Hq). For each n e N there exists a linear projection 
operator Hq : L^(il) Q" such that, for ah s' e (l,co), Hq is stable. In other 
words, there exists a constant c^' > 0, independent of n, such that 

The stability of Hq and the approximation properties of Q" c L" (f2), stated in 
Assumption [5j imply that Hq satisfies 

||g — nQqjls' ^0, as ri ^ CX) for all q e (fi) and s' e (1, oo). (3.5) 

As a consequence of (|2.ip and Assumption [5] (compare also with (|3.3p ) the fol- 
lowing discrete counterpart of (j2.ip holds; see, |BBDR10| . 

Proposition 8 (Inf-sup stability). For all s, s' e (1, co) with 7 + p- = 1, there exists 
a constant /3s > 0, independent of n, such that 

sup ^'^^Ly!,' '^^^ > I3s II g lis' for all Q e QJJ and all neN. 

0#VeV" ll»'l|l,s 

This property is in turn equivalent to the existence of a discrete Bogovskii oper- 
ator; compare with |BF91[ §11, Proposition 1.2]. 



8 



L. DIENING, C. KREUZER, AND E. SULI 



Corollary 9 (Discrete Bogovskii operator). Under the conditions of this section, 
there exists a linear operator : divV" — > V" and a positive constant /3g, inde- 
pendent of n, such that 



div(<B"i?) = H and ||<B"ff||i.s ^ sup 



WQWs' 

for all He divV", n e N. 

It follows from Corollary E] by Holder's inequality that /3^. ||«8"i7|| < \\H\\s. 
However, we shall need in the proof of Lemma [50] the stronger statement from 
Corollary [9] 

3.2. Discretely divergence- free finite elements. As in |Tem84) we wish to en- 
sure that the discrete counterpart of the convection term inherits the skew-symmetry 
of the convection term. In particular, upon integration by parts, it follows that 

- I {v®w) -.Vhdx = I {v®h) -.Vw + {d\vv){w ■h)dx (3.6) 

for all v,w,h e P(O)'^. The last term vanishes provided that divv = 0, and then 

/ {v ®v) : Vv dx = 0. 
Jn 

It can be easily seen that this is not generally true for finite element functions 
V", even if 

<divy, Q^ = forallQeQ", 
i.e., if V is discretely divergence-free. However, we observe from (|3.6p that 

- / {v(g)w) : Whdx =\ / {v®h) -.Vw - {v®w) : Whdx =: B[v, w, h] (3.7) 
Jn Jn 

for all v,w,he W^'^^^{n)'^. We extend this definition to W^-'^^ilY in the obvious 
way and deduce that 

B[v, v,v] = for all v e W^'^^inf. (3.8) 

We further investigate this modified convection term for fixed r,r' e (1, '») with 
i + p- = 1; recall the definition of r from (|l.ip . We note that f > 1 is equivalent to 
the condition r > In this case we can define its dual f' 6 (1, cx)) by i + J7 = 1 
and we note that the Sobolev embedding 

W^-^'iny ^ L^^i^Y (3.9) 

holds. This is a crucial condition in the continuous problem, which guarantees 

/ {v®w) -.Vhdx ^ c\v\i^r\w\i,r\h\i,f' (3.10) 
Jn 

for all v,w,h(i W'^^''='{nY] see |BGMSG09| and Section below. Because of the 
extension p.7p of the convection term to functions that are not necessarily point- 
wise divergence-free, we have to adopt the following stronger condition in order 
to ensure that the trilinear form is bounded on W^'^{^Y ^ W^'^{^lY ^ 

W^'^ {^Y- In particular, let r > in order to ensure that there exists s e (1, co) 
such that 7 + ^ + 7 = 1- In other words, we have for v,w,h e W'^''^{ftY that 

/ (divt;) {w ■ h) dx || div 7;||r||ii'||2f < c ||v||i^r||'U'||i,r||/j'||i,f', 
Jn 

with a constant c depending on r, VL and d. Here we have used the embeddings 
and Wl'^\nY ^ L^iyty. Consequently, together with ([XTU)) we thus obtain 

B[V, W, h] ^ c||tl||i,r||l«||l,r||^||l,f'- (3.11) 
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Example 10. In [BBDRIO] it is shown that Assumptions [6] and [7] are satisfied by 
the following velocity-pressure pairs of finite elements: 

• The conforming Crouzeix-Raviart Stokes element, i.e., continuous piecewise 
quadratic plus bubble velocity and discontinuous piecewise linear pressure ap- 
proximations (compare e.g. with |BF91[ §VI Example 3.6]); 

• The Mini element; see, |BF91[ §VI Example 3.7]; 

• The spaces of continuous piecewise quadratic elements for the velocity and piece- 
wise constants for the pressure f |BF91[ §VI Example 3.6]); 

Moreover, it is stated without proof in [BBDRlO] that the lowest order Taylor-Hood 
element also satisfies Assumptions [S] and [T] 

3.3. Exactly divergence-free finite elements. Another way of retaining the 
skew-symmetry of the convection term and ensuring that p.8p holds is to use an 
exactly divergence-free finite element approximation of the velocity. In addition to 
Assumptions ini and [7] in Section l3.1l we suppose that the following condition holds. 

Assumption 11. The finite element spaces defined in Section I57T] satisfv 

divV"czQf,\ neE. 

This inclusion obviously implies that discretely divergence-free functions are auto- 
matically exactly divergence-free, i. e., 

V^i^ = {V eY'': divV = 0}, ne N. 
According to (|3.6p . in this case, we define 

B[v,w,h]:=- {v0w):\'hdx (3.12) 
Jn 

for aU v,w,he W^'^^iflY and obtain 

B[v, v,v] = for all v e W^;^^^{nf. (3.13) 

Recalling (|3.10p , with Assumption [TI] the convection term can be controlled under 
the weaker restriction r > i. e., for v,w,h e W^'^'{^)'^ ^ have that 

B[V, W, h] ^c\\v\\i,rMi,r\\h\\l,r', (3.14) 

where, as before, i + = 1 with f from p.ip . The constant c > only depends 
on r, n and d. 

Admittedly, finite element spaces that simultaneously satisfy Assumptions |6l [7] 
and [TT] are not very common. Most constructions of exactly divergence-free finite 
element spaces in the literature are not very practical in that they require a suf- 
ficiently high polynomial degree and/or restrictions on the geometry of the mesh; 
see |AQ92[ ISV85[ |QZ07[ IZhaOS] . In a very recent work [GNllj Guzman and Neilan 
proposed inf-sup stable finite element pairs in two space-dimensions, which admit 
exactly divergence-free velocity approximations for r = 2. A generalization of the 
Guzman-Neilan elements to three dimensions is in progress |GN12| . We shall show 
below that the lowest order spaces introduced in |GN11| simultaneously satisfy As- 
sumptions [HI [7] and [TT] for r 6 [ 1 , cx) ) . 

Example 12 (Guzman-Neilan elements |GN11| ). We consider the finite element 
spaces introduced by Guzman and Neilan in |GN11[ Section 3] on simplicial trian- 
gulations of a bounded open polygonal domain f2 in R^. In particular, we define 

Pv := Fi(i?) © span{curl(&j) : i = 1,2,3} © span{curl(i3j) : i = 1,2,3}. 

Here Fi{E) denotes the space of affine vector- valued functions over E. Let further 
{Aj;}i=i^2,3 be the barycentric coordinates on E associated with the three vertices 
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{zi}i=i,2,3 of E, i.e., Xi{zj) = Sij. Then, for i e {1,2,3}, we set bi := and 
Bj denotes the rational bubble function 



B, 



X X2 \2 



(A, + Aj+i)(A, + \,+2) 

which can be continuously extended by zero at z^+i and Zi+2; the index i has to be 
understood modulo 3. Thanks to properties of the curl operator, the local pressure 
space 

Pq := divPv 

is the space of constant functions over E. 

It is clear from |GN11| that the related pairs of spaces {V",Q"}„eN (compare 
with (|3.ip ) satisfy Assumption [TT] For IIq we can use a Clement type interpolation 
or simply the best-approximation in Q" ; clearly, both satisfy Assumption [71 The 
approximability assumption. Assumption [51 follows with the mesh-size tending to 
zero. It remains to verify Assumption [51 To this end we analyze the interpolation 
operator proposed in |GN11| . In particular, let Eg : WQ'^{n)'^ L" be the Scott- 
Zhang interpolant |SZ90j into the linear Lagrange finite element space L" := h{Qn) 
over the shape-regular triangulation Qn of il, n e N. Then II^j^ : VFq' (il)'^ V" 
is defined by 

(nSi.t;)(z) := (n»(z), ^eAA", 



Is Js 

where A/"" and 5" denote the vertices, respectively edges, of the triangulation On, 
n 6 N. This operator is a projector and thanks to [GNlli (3.14)] and the fact that 
L" c V" it thus remains to prove the stability estimate (|3.2p in Assumption [51 

To this end we fix n e N. Although the claim can be proved using the techniques 
in |GN11| , this would necessitate the introduction of additional notation. Thus, for 
the sake of brevity of the presentation, we give an alternative proof. According to 
[GN11| the interpolation operator II^j^ is correctly defined by p.lSp . Let {zi}i=i.2.3 
be the vertices of E and let {S'i}i=i_2,3 be its edges. Then, any function V e Py is 
uniquely defined by V{zi) and /a. Vds, i = 1,2,3. This implies that the mapping 



3 „ 

V|y(z,)|+ / Vds 

,=1 Js, 



Ve 



is a norm on Py. Hence, equivalence of norms on finite-dimensional spaces together 
with (|XT5|) yield that 

3 

/ \U2;^v oF-^'\dx^cj] (\n2,^v o F-/{z,)\ + I / n^.^v o F'/ ds\) 

JE J Si 

3 „ 

= c;^ (|n>oF^i(2,)| + | / voFeUs\), vew'^-\nf. 
i=l J Si 

By a scaled trace theorem and properties of the Scott-Zhang operator we arrive at 

[ \S^^^v\ dx ^ c [ \v\ + \E\^\Vv\dx, veWl'\^f. 
J E inj 

Note that : Wl'^{nY ^ V" is a projector and that L" c V". Thus the 

inequality p.2p follows from a standard inverse estimate and the Bramble-Hilbert 
Lemma; compare also with [BBDRlOl Theorem 3.5]. 
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3.4. The Galerkin approximation. We are now ready to state the discrete prob- 
lem. Let {V",Q"}„gN be the finite element spaces of Section [3.21 or 15751 and let 
B : W^^'^iny X Wg^'^inY x W^^'^inY ^ R be defined correspondingly. 

Forne Nwecallatripleof functions ([/", P", S"{DU")) 6 V"xQgx L'^'(rj) 
a Galerkin approximation of p.2p if it satisfies 

/ S"{-,DU'') : DVdx + B[U'\ U'\ V] -(divV, P"^ = </, V) 
Jq 

for all Ve V", 

[ Qdivi7"da; = for all Q e Q". 
Jn 

(3.16) 

Restricting the test-functions to V^j^ the discrete problem p.l6p reduces to find- 
ing i7" e V^i^ such that 

[ S'"(-, DU") ■.DVdx + B[U", [/", V] = </, Vy for ah V e V^^. (3.17) 

Thanks to p.Sp . respectively (|3.13p . it follows from Lemma [T] and Korn's inequal- 
ity (|2.2p that the nonlinear operator defined by the left-hand side of p.l7p is coercive 
on V^jy. Since the dimension of VJJ;^ is finite, Brouwer's fixed point theorem ensures 
the existence of a solution to (j3.17p . The existence of a solution triple to p.l6p then 
follows by the discrete inf-sup stability implied by Proposition [8] 

Of course, because of the weak assumptions in the definition of the maximal 
monotone r-graph, p.l6p does not define the Galerkin approximation i7" uniquely. 
However for each n e N we may select an arbitrary one among possibly infinitely 
many solution triples and thus obtain a sequence 

{(C/",P",5"(.,DJ7"))}„^^. (3.18) 

3.5. Discrete Lipschitz truncation. In this section we shall present a modified 
Lipschitz truncation, which acts on finite element spaces. This discrete Lipschitz 
truncation is basically a composition of a continuous Lipschitz truncation and the 
projector from Assumption [51 For this reason we first introduce a new Lipschitz 
truncation on Wq'^{VlY, based on the results in |DMS08I IBDF12I IBDS12) . which 
provides finer estimates than the original Lipschitz truncation technique proposed 
by Acerbi and Fusco in |AF88) . 

For V e L^{W^) we define the Hardy-Littlewood maximal function 

M{v){x) := sup f \v\ dy. (3.19) 

For s e (1, oo] the Hardy-Littlewood maximal operator M is continuous from L'^{M.'^) 
to L'^{M.'^), i. e., there exists a constant Cg > such that 

||M(i;)|U.(Kd) s; cs Ml^^m") for aU v e L'iR''), (3.20) 

and it is of weak type (1, 1), i. e., there exists a constant ci > such that 

sup A \{M{v) > X}\ ci||?;||ii(Rd) for aU v e L^R'^); (3.21) 

A>0 

see, e.g., ^0^. For any v e W^-\R'^Y we set M{v) := M(|t!|) and M{Vv) := 
M(|Vi;|). 

Let V 6 VFq'^(51)'^; we may then assume that v e W^'^(W^Y by extending v by 
zero outside il. For fixed A > we define 



Ux{v) := {M{Vv) > A} , 



(3.22a) 
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and 

nx{v) := R\(Uxiv) nn) = {M{Vv) A} u (3.22b) 

Since M{\7v) is lower semi-continuous, the set U\{v) is open and the set T-L\{v) 
is closed. According to |DMS08) it follows that v restricted to 'H\{v) is Lipschitz 
continuous and therefore also bounded. More precisely, we have that 

\v{x)-v{y)\^c\\x-y\ (3.23) 

for all x,y e TL\{v), where the constant c > depends on ft. 

It remains to extend ^i-^^^^,-) to a Lipschitz continuous function on M.'^. The result 
in [DMS08) is based on the so-called Kirszbraun extension theorem (cf. Theorem 
2.10.43 in |Fed69| ) and uses an additional truncation of v with respect to M{v). 
This can be avoided by proceeding similarly as in |BDF12[ IBDS12| . i. e., extending 
'^I'Hxiv) by means of a partition of unity on a Whitney covering of the open and 
bounded set U\(v). To this end, we assume w.l.o.g. that U\{v) ^ 0; otherwise v 
does not need to be extended since 'H\{v) = K''. According to |G04i |BDF12| there 
exists a decomposition of the open set U\(v) into a family of (dyadic) closed cubes 
{QjljEN, with side lengths Ij := ^(Qj), j 6 N, such that 
(Wl) UjeN Q-j ~ ^a(i') and the Qj's have pair-wise disjoint interiors. 

(W2) 8Vd£(Qj) d\st{Qj,eUx{v)) ^ 32Vd£{Qj). 
(W3) If Qj nQk^ for some j, keN, then 

<g <g 2. 



2 

id _ 1 \od 



(W4) For a given Qj there exists at most (S^* - 1)2'* cubes Qk with Qj n Qk ¥= 
For a fixed cube Q e R*^ with barycenter z and any c > 0, we define 



cQ:=ixeW: max \x^ - zA ^ ce(Q) 

[ 1=1, ...,d 

Hence, it follows from |(W2)| with 6*^ := 2 + 64\/d, that 

(OdQj) r^ Hx{v) ^ for all j e N. (3.24) 

Let Q* := ^J^Qj and Q** := fQj. Thanks to pV4)| the enlarged cubes Q**, 
j 6 N, are locally finite, i.e.. they satisfy Y! ■ Xo** ^ ^ with a constant c > 0, which 
depends on d only. Thanks to the overlaps of the Q*'s, there exists a partition of 
unity {T/ijIjeN subordinated to the family {Q*}jeN with the following properties: 

• Yjj V'j = XUx(v) and ^/'^ 1 for all j e N; 

• X|c, < V'j =S Xq* , for all j e N; 

• e Co°°(Q*) with IW'jl ^cij\ for all j e N. 

The Lipschitz truncation of v is then denoted by vx and is defined by 

^^^^fS.eM'A,^, n^l^xiv), ^3 25a) 

) elsewhere, 

where 

vdx ifQ*crj, 

^ (3.25b) 
elsewhere. 

We emphasize that the definition of the functions Vj, j e N, here differs from 
the one in |BDF12| . since we need to preserve the no-slip boundary condition for 
the velocity field on under Lipschitz truncation. Combining the techniques 
of |DMS08) and |BDF12| we obtain the following result; for ease of readability of 
the main body of the paper, the proof of Theorem is deferred to the Appendix. 
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Theorem 13. Let A > and v e VFq' (ri)"*. Then, the Lipschitz truncation defined 
in p.25p has the following properties: V\ e W^^'^'iri)'^ , and 

(a) vx = V on Hx{v), i. e., {v ^ vx} <= Ux{v) r^ n = {M{\7v) > A} n r2; 

(b) II^aIIs c \\v\\s for all s e [1, co], with v e L'^{Q,y; 

(c) ||Vi;a||,s cIlV'ull,, for all s e [l,co], with v e W^'''{nY; 

(d) |Vva| 5S cXxux(v)nn + |Vf | X-HxM almost everywhere in M."^. 

The constants appearing in the inequalities stated in parts (b)[ (c) and |(d)| depend 
on Q and d. /w |(b)| and (c) they additionally depend on s. 

We next modify the Lipschitz truncation so that for finite element functions the 
truncation is again a finite element function. To this end we recall the definition of 
the finite element space V" = V(C/„) of Section [XT] or [XHl 

Let A > and fix n e N. Since V" c; W^'^^n)'^, we could apply the truncation 
defined in p.25p . However, since in general the Lipschitz truncation Va of V e V" 
does not belong to V", we shall define the discrete Lipschitz truncation by 

y„,A := n^iv^A e V". (3.26) 

According to the next lemma the interpolation operator H^i^ is local, in the sense 
that it modifies Vx in a neighborhood of Ux{V) only. 

Lemma 14. Let V e V". With the notation adopted in this section, we have that 

{V„,x ^V}^ ni{V) := interior (|J {n]^ \ E e g„ wtth E n Ux{V) ^ (2 

Proof. The stated inclusion follows from p.2p in Assumption [51 In particular, 
let E e Gn he such that E c K'^\r25;(y); then, n% c TixiV). Consequently, 
by Theorem [T^aJI we have Vx = V on fi^. Hence we deduce from p.2p . our 
assumption that V e ¥", and the fact that H^j.^, is a projector, that 

-f \V- n^wVxl dx = / |HSi,(F -Vx)\dx^C'f \V- Vx\ dx = 0, 
Je Je Jn^ 

i.e., V = Vn,x = Ildiv^-^ ^- This proves the assertion. □ 

The set rt^{V) from Lemma [TTl is larger than Ux{V) n il. However, the next 
result states that we can keep the increase of the set under control. This is the key 
observation for the construction of the discrete Lipschitz truncation. 

Lemma 15. For n e N, V e ¥" and A > 0, let f2^(V) be defined as in Lemma\T^ 
Then, there exists a k e (0, 1) only depending on Py and the shape-regularity of Qn, 
such that 

Ux{V) n O c: ni{V) cz U^x{V) n n. 

Proof. Thanks to the definition of fi^KV), the first inclusion is clear. It thus remains 
to show the second inclusion. In order to avoid problems at the boundary dfl we 
extend V to R'' by zero outside $7. Let x e rt^{V); then, there exists E e Qn, 
E n Ux{V) such that x e fl"^. Consequently, by p.22ap and ([^1^)) . there 
exists an xq e and an i? > such that 

/ |VF| dy > A. 

JBr{xo) 

Suppose that Bj^fxo) c (Jl^ u (M''\n)); then, thanks to norm-equivalence in 
finite-dimensional spaces, we have, by a standard scaling argument, that 



A < 



/ |VF| dy^ ||Vy|Ux(o^) ^Si/ |vy|dy, 
Jbr{xo) Jn^ 
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where the constant ci depends solely on Py and the shape- regularity of Q. Let 
Bp{x) be the smallest ball with Bp{x) zd and observe that |i3p,(a;)| ^ C2 
with a constant C2 > depending only on the shape- regularity of Q. Consequently, 

Af(Vy)(x) Ss / |Vy| dy^^i \VV\ dy ^ -J-^X. 

Jb,{z) C2 Jit'^ C1C2 

In other words, we have that x e Ui^£-^co)-^x{^) ^■ 

We now consider the case Bb{x{)) <^ (jl^ u (R''\57)). Since xa e E, it follows 
that C3R 5= diam(£'), with a constant C3 > only depending on the shape- regularity 
ofQn- As a; e fi^. there exists a constant C4 > 1 such that Bc^j^{x) zd {Eyj Bji{xQ)). 
Hence, 

M{VV){x)^'f |VV| dy £4 ^ / |VV| dy > £4 ^A, 

and we deduce that x e U^-d^{V) n il. Combining the two cases, the claim follows 
with K := min{(ciC2)^^, £4^''}. □ 

We are now ready to state the following analogue of Theorem [T^] for the discrete 
Lipschitz truncation p.26p . 

Theorem 16. Let A > 0, n e N and V e V". Then, the Lipschitz truncation 
defined in p.26p satisfies Vn,\ e V", and the following statements hold: 

(a) F„,A = V onR''\n''^{V); 

(b) \\Vn,x\\i,s iic\\V\\i,s /or 1< s CO; 

(c) |VF„,a| s; cAxnj(v) + |Vy| XEd\j2j(v) < cA almost everywhere in R''. 



The constants c appearing in the inequalities in parts [(b)] and (c) depend on fl, d, 
Py and the shape-regularity of {Qn\n€n- fa |(b)| the constant c also depends on s. 



Proof. Assertion (a) is proved in Lemma (TSl Estimate |(b)| is a consequence of 
Theorem [T51 and the ly^'^-stability of nji^; compare with (|3.5|) . The bound (c) 



follows from Theorem [I^d)] and the W^'*' stability of H^j^; see (031). □ 

The following corollary is an application of the discrete Lipschitz truncation to 
(weak) null sequences. It is similar to the results in [DMS08J and jBDF12| . Its 
analogue for Sobolev functions is stated in Corollary in the Appendix. 

Corollary 17. Let 1 < s < co and let {E'^}nm c wI''{VlY be a sequence, which 
converges to zero weakly in Wg^'^(f2)'', as ?i — > cx). 

Then, there exists a sequence {A„j }„j(=n c ^ with 2^^ ^ A„j ^ 2^^^ such that 
the Lipschitz truncations E^^'-' := a„ j ; J ^ have the following properties: 

(a) E""'^ 6 V" and E""'^ = E"" on R'^\ni^ ^{E"") ; 

(b) ||V£;"'^||oo<cA„,,; 

(c) ^ m L'^in)"^ asn^ ao; 

(d) VE""^ — * in L'^inY'"^ asn^ co; 

(e) For all n,j e N we have ||A„,-,xoj ,{E")h < c2-^ || V£;"||s. 



The constants c appearing in the inequalities \{h)\ and (e) depend on d, Q, Py and 
the shape-regularity of {Qn}nen- The constant in part (e) also depends on s. 
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Proof. We first construct the sequence A„ j- and prove (e) Let k > be the constant 
in Lemma [ini Then, for g e L'*(M''), we have 



[ \g\'dx=[ r\'sf-\^gi^^t^dtdx^ [ Vn'sr 

JR-' JR-^ Jo JR'' 



X{|gl>A.2" + i} ' 



2^ + 1-1 

^J] J] A.-^s2™-^|{|g|>2™+i}|- 

We apply this estimate to g = MiVE^) and use the boundedness of the maximal 
operator M (cf. (j3.20p ) to obtain 

2^ + 1-1 

k'Yj Yj |{M(V£;") > 2™+i}| < c, \\M{VE")\\l < c,\\VE"\\l. 

For fixed n, j the sum over m involves 2-' summands. Consequently, there exists an 
integer A„ j e {2^^ , . . . , 2^^ such that 

A^,^^. \{M{VE") > kA„,,}| < 2-^n-'cs\\VE"rs. 

This, together with Lemma [TSl proves (e) Assertions (a) a nd |(b)| are then direct 
consequences of Theorem [T^a)| and (b) It remains to prove (c) and |(d)[ 



Let j e N be fixed. From [(b)] we deduce that the sequence is bounded 

in Wq''^ {fl)'^ . Hence, there exists a (not relabeled) subsequence such that V.E"'-' 
converges weak-* in L^'(f2)'*'*'^ and E"'-' converges strongly in L''^{il)'^ as n ^ cx). 
Thus we have proved part (c). 

To prove (d), we proceed as follows. Thanks to the uniqueness of the limits, it 
suffices to prove that E'^'-' weakly in Wf^'^{ftY. To this end, we note that the 
compact embedding Wq''^{Q,)'^ L*(51)'' implies that 

E" in L^iny as n ^ cx). 
Let {£^A„ j }neN be the sequence of Lipschitz-truncated functions, defined according 
to ((X^ . Then, thanks to the boundedness of {£;"}„eN in VKJ '"(r?)'', TheoremH^c)] 
and |(b)[ we have that 

El^ , ^ weakly in W^'^n/, as n ^ oo. 

Thanks to the definition of the discrete Lipschitz truncation in p.26p , the desired 
assertion follows by observing that the norm-continuity (|3.3p of HJJj^ in Wo'*(ri)'' is 
equivalent to continuity of njj;^ in the weak topology of WQ''^{fl)'^ (see, for example. 
Theorem 6.17 in |AB2006| ). □ 

4. The main theorem 

After the preceding considerations, we are now ready to state our main result. 
Its proof is presented in subsections ij4.1L t i4.41 

Theorem 18. Let {V",Q"}„en be the sequence of finite element space pairs from 
Section\KE (respectively and let {{U" , P" , S"{-, DU"))}^^^^ be the sequence 
of discrete solution triples to p.l6|) constructed in p.lSp . 

If r > (respectively r > -j^), then there exists a solution {u,p,S) e 

Wo'^'iny X Lf^ifl) X L'''{ny'"^ of ((L2l), such that, for a (not relabeled) subsequence, 
we have 

U" — u weakly in W^'^'iflf, 
P" ^ p weakly in Llf(^l), 
S"i-,DU'') - S weakly m U' {^f'^'^ . 
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4.1. Convergence of the finite element approximations. We begin the proof 
of Theorem [TS] by showing the existence of a weak hmit for the sequence of solution 
triples. 

Lemma 19. Let {V",Q"}„eN be the sequence of finite element space pairs from 
Section [KE (respectively [j;^ and let {([/", P", £)[/")) be the sequence 
of discrete solution triples to (j3.16p constructed in (|3.18p . 

If r > (respectively r > ^^), then there exists {u,p,S) e W^'^^il)'^ x 

L5(f2) X (17)''^'', such that, for a (not relabeled) subsequence, we have 

t/" — M weakly in W^'''{nf, 

P" ^ p weakly in L^^fl), 

S"{-,DU'') — S weakly in L''' {VL)'^'"^ . 

Moreover, the triple {u,p,S) satisfies 

I S : Dv- {u(^u) -.Vvdix-idivv, p)^^ = (f, v), for all v e Wq'^' {n^ 
'''' ^ (4.1) 

/ g div M da; = 0, for all q e U {Vl) . 

Proof. We divide the proof into four steps. 

Step 1: From (|3.16p we see that i7" is discretely divergence- free and thus, thanks 
to (|3.17p and p.Sp (respectively (|3.13p ). we have that 



/ 5"(-,D(7") :!?[/" dx = </, i7">< ||/||_i,,,I|[7"| 
Jn 



l,r- 



The coercivity of 5" (Lemma[T]) and Korn's inequality (|2.2p imply that the sequence 
{U"}nen c Wq'^(P,Y is bounded, independent of n e N. This in turn implies, again 
by Lemma [U the boundedness of {5'" (DL''")}„(=n hi L''' (D.)'^^'^. In other words, 
there exists a constant c > such that 

||Lr"i|i,,. + ||5'"(-,DLr")||,., c, for all n e N. (4.2) 

As r 6 (l,co), the spaces Wq^'' (fi)'' and {VtY^'^ are reflexive and thus for a (not 
relabeled) subsequence there exist u e W^'^ {Vl)'^ and S e U" {VlY^'^, such that 

Lf" — M weakly in Wo''{9)'^ (4.3) 

and 

S"{-,DU") — S weakly in L''' {n)'^'''\ (4.4) 
as n ^ <». Moreover, using compact enibeddings of Sobolev spaces, we have that 



V ' d-rj ' ' (4.5) 

s e (l,oo), otherwise. 

Thanks to (j4.3p we have by p.5p . for arbitrary q e L'^ (fl), that 

0=1 {U^q)divU"dx ^ [ gdivitdx, (4.6) 
Jn Jn 

i.e., the function u e Wq'^ (fl)'^ is exactly divergence-free. 

Step 2: Next, we investigate the convection term. Let v e Wq''^(JI)''' be arbi- 
trary and define V" := ITliyV. We show that 

B[U'\ U'\ V"] [ {u®u): Vi;dx. (4.7) 

Jn 
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Thanks to the assumption r > and (j4.5p . it foUows that 

J7" (X) [/" ^ M (X) M in L-'inY""^ for aU s6 [l,f), 

with f > 1 as in ([LT|) . By (j^ . we have that F" ^ in W^^ s' e (f',cx)), 
and hence we obtain that, as n — > co, 



(4.8) 



- / ((7" (X) (7") : VV" dx ^ - /" (m (X) It) ; dx. 
Jn Jn 

This proves (j4.7p for the exactly divergence-free approximations from Section 13.31 
We emphasize that we have only required so far that r > 

In order to prove (|4.7p for the finite element spaces of Section 13.21 and thus for 
the modified convection term defined in p.7p . we recall from p.6p that 

f ([/" ® F") : VU" dx = - f (i7" g) [/") : VV" + (div [/") (7" • V" dx. 
in Jn 

For the first term we have already shown convergence in (j4.8p . In view of the 
definition of B in (|3.7p it thus remains to prove that the second term vanishes in 
the limit n ^ co. To this end, we observe by (j4.5p and Assumption [6] that 



strongly in for all 



se(^l,^j, ifr<d, 
,se(l,cx)), otherwise. 



Thanks to the stronger restriction r > now, this last statement holds in par- 
ticular for s = r' . Hence, together with (|4.3p and (|4.6p . we deduce that 



[ (div (7") i7" • dx ^ /" (div m) m • dx = 
Jo. Jn 



as n ^ cx). 

Step 3: We combine the above results. Recall that by (j4.6p we have divM = 0, 
which is the second equation in (j4.ip . For an arbitrary v e WQ div(^)'^ 
njJi^v, n e N. Thanks to ([231), we have that F" e VJf;^ and V" ^ t; in l^o'''(17)'^ 
for aU s e (l,cx)). Therefore, using (|4.3p . (|4.4p and (|4.7p . we obtain 

/" 5"(-, r>J7") : DV" dx + B[U''\ U'\ V"] = </, V") 
i i 

/ S : Dv + div(M ®u) ■ vdx (/, v) 

Jn 

as n ^ (X>. Since S e L'''(rj)''^'', / e W-^''''{ny and u®u e L^(17)'^^'^, by a 
density argument, we arrive at 

/ S : Dv + div{u 0u) ■ vdx = (/, v} 
Jn 

for a\lveW^:^,,{nr. 

Step 4: We now prove convergence of the pressure. Thanks to the restriction 
r > we have, as in p. lip , that 

<divF, P"^ = / S'\-,DU") : DVdx + B[U'\ U" , V]-(f, V} 
Jn 

s; \\S^{-,DU")\UDV\\r + c\\UXjV\U^f' + \\f\\-iAV\\i,r 

for all V 6 V". By (|4.2p and the discrete inf-sup condition stated in Proposition[51 it 
follows that the sequence {P"}„(=n is bounded in the reflexive Banach space io(f^). 
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Hence, there exists p e ^^(f^) such that, for a (not relabeled) subsequence, P" ^ p 
weakly in L^iyi). On the other hand we deduce for an arbitrary v e Wq '°°(51)'' that 

<div^, P">^ = (divnSi,.;, P">,, + <div(t; - nSi,^), P"X, 
S : Dv + div{u ®u) -vdx- </, v} + 

Jn 

as n — > 00, where we have used (|4.4p . (|4.7p . p.4p and the boundedness of the 
sequence {P^jneN in io(^)- This completes the proof of the lemma. □ 

For the main result. Theorem I18[ it remains to prove that 

{Du{x), S{x)) e A{x) (4.9) 

for almost every x e il. The proof of this is the subject of the rest of Section 2] 

4.2. Identification of the limits. In this section we shall first briefly discuss 
properties of the maximal monotone r-graph introduced in (Al) (A5) Here we 
follow the presentation in p3GMSG09j . Application of the fundamental theorem on 
Young measures (cf. Theorem [3]) leads to a representation of weak limits, which is 
a crucial step in proving (|4.9p . 

According to |FMT04| there exists a function : 51 x R^ym such that 

A{x) = {{d, ct) e M^-^ X K^;^^ -.a d = <t>{x, a + 5)} , (4.10) 

and 

(a) 0(a;, 0) = for almost every x e fi; 

(b) (/'(•, x) is measurable for all x ^ ^iyit'^ 

(c) for almost all x e the mapping 4>{x, ■) is 1-Lipschitz continuous; 

(d) the functions s, d : fl x W^^^^ Rfy'<„f , defined as 

six,x)-=Hx + <l)ix,x)), d{x,x) ■■= ^{X-Hx,x)) (4.11) 

^dxd 
sym ' 



satisfy, for almost every x e and all x 6 the estimate 



s{x, x) ■■ d{x, x) ^ -mix) + c [\d{x, xW + \s{x, x)\ 

We emphasize that this is in fact a characterization of maximal monotone r-graphs 
A satisfying (Al) (A5) without the second part of |(A2)[ 

We recall the selection from Section and, as in |BGMSG09) . we define 

bn{x) := [ {S*{x, C) - S*{x, Du)) : (C - Du{x)) dy.^^{C). (4.12) 

where we have used the abbreviation /i" := /^^^j-^.). The next result, whose proof 
is postponed to the next section, states that 6„ vanishes in measure. 

Lemma 20. With the definitions of this section we have that 6„ Q in measure. 

Actually, employing the above characterization of A^ the limit of the sequence 
{6„}„gN can be identified in another way by using the theory of Young measures. 
To this end we introduce 

G,(C) := S*{x, C) + C, xen,C^ <m, (4.13) 
and define the push-forward measure of the measure /i^ from (j2.5p by setting 

<c(C) = Mc(G^'(C)) for aU C e ^{Ky^). (4.14) 
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We recall from ^2.^l^al) that S* is measurable with respect to the product cr- 



algebra 2{n) ® «8(R^^, and therefore so is Gx- Consequently, the measure i/"> 
is well-defined and, thanks to properties of the mollifier 77" from Section 12.21 it is a 
probability measure. From the definitions of the functions s and d it follows that 
S*{x, C) = s{x, Gx{C)) and C = d{x, Gx{C))- We thus have, by simple substitution, 
the identities 

S-{x,DU'^)=[ s(x,C)dz.:(C), (4.15a) 
DU^{x) = [ d{xX)AK{0, (4.15b) 

as well as 

bn{x)= I {s{xX)- S''{x,Du{x))) ■.{d{xX)- Du{x))dv'^{C), (4.15c) 

where we have used the abbreviation i^" := ]^■^n^^x)^ 

In order to identify the limit we apply the generalized version of the classical fun- 
damental theorem on Young measures stated in Theorem |31 Recall from Section [^751 
that L'^{n-m{Ri^,^)) is a separable Banach space with predual i^^^, C'o(Rfym))- 
For every 71 e N the mapping x ^ belongs to L'^{^] 9Jt(Rfy„f)). To see this let 
g 6 L^{Vl; Co(Msym)). In view of (|4J4)) and ((2?5|) . a simple substitution yields 

/ g{x, C) du^X) = I v^{DU'\x) - C) g(x, G,(C)) dC. 

It remains to prove that x ^ /u^xd rpiDV^ix) — C) .9(2;, Gx{C)) dC is measurable 
and integrable. It follows from (A5) the definition (|4.13p of Gx and property |(al)| of 



S* that h is £(fi)(x)*8(Rfjf„f ) measurable. Moreover, 77", 71 5= 1, are smooth functions 
and g e L^(fi, Co(Rsy^nj)), and therefore integrability follows from Fubini's theorem. 

Thanks to the properties of 77" it is clear that is a probability measure a.e. 
in 51. Hence ||^'"|lia:(Q.<j)^(R<i>''')) = 1 and thus the sequence {j^'^jneN is bounded 
in L2; (0; 93t(Ms^j^)). Therefore, by the Banach-Alaoglu theorem, there exists v e 
L^(5I; 9Jl(Rfy^j^)) such that, for a (not relabeled) subsequence, 

v'' ^* V weak-* in Li^'(57;SH(R^y''^)). (4.16) 

The next Lemma follows from the generalization of the fundamental theorem 
on Young measures from |Gwi05| (see Theorem and Chacon's biting lemma 
(Lemma 0]). Its proof is postponed to Section 

Lemma 21. With the notations of this section, Vx is a probability measure a.e. 
in n and there exists a nonincreasing sequence of measurable subsets Ek c 51, with 
\Ek \ 0, such that for all k e N we have for a (not relabeled) subsequence that 

bn{x) - [ {s{x, C) - S*{x, 0) : {d{x, C) - Du{x)) dux{C) =: b{x) (4.17) 

weakly in L^{Q,\Ek) as n ^ co. Moreover, we have that 

S{x) = / s{xX)dvx{C) and Du{x) = / d{xX)di^x{C)- 



sym ^^sym. 



We deduce from (|4.9p and Lemma [5T] that, to complete the proof of Theorem [T51 
we need to show that 

[Duix), s{xX)dux{C)] e A{x) for a.e. x e fl. (4.18) 
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This follows from the two preceding Lemmas exactly as in |BGMSG09[ p.l31fF]. 
To be more precise, the proof is based on noting that for each S e R''^'' the set 

Cg := {(T 6 W^^^ : {5, a) e A{x)} is convex for a.e. x e n; (4.19) 

compare with |(A2)| and |(A3)[ Combining Lemmas [50] and [5T] we deduce that b = 
a.e. in fl. Hence it follows from |(A2)] that 

{C6Kf;,„^: {s{x,C)-S''{x,Du{x)) : {d{x,C) - Du{x)) > 0} suppi.,. 

We split supp Ux into the two sets 

uji{x) := e suppi/x- s{x,<^) = S* {x, Du{x))} and uj2{x) := supp z^^r \ wi (x) . 

We investigate the pairing in (|4.18p on the two sets uji{x) and uj2{x) separately. On 
uj2{x) we have by |(A2)| that d{x,C) = Du{x). Therefore, on noting that ^ '(^^^i)) 
is a probability measure on ijJ2 {x) , one can show with (|4.10p , (|4.1ip and (|4.19p that 



[du{x), f s(x,C)cl(- 



e A{x) for a.e. x eVL. 

On the other hand it follows from the definition of cji [x) that 

s{xX)<^Vx{C) = Vx{^i(:x))S*{x,Du{x)). 

Thanks to the properties of S'*, we have that (^Du{x), S* {x, Du{x))) e A{x); 
compare with Section 12.21 Using the fact that v^. is a probability measure, we 
deduce that 

/ s{xX)dMC)= [ s{xX)dMC)+ [ s{x,C)di^x{C) 

■/R?ym •'^1(2;) JUJ2{X) 

= ,^xiu;,ix))S*ix,Duix)) + MM^)) f <^X)d{^;^^^^) 

J u:2(x) 

is a convex combination of functions. Moreover, due to the above observations, for 
a.e. X e Q, each of the two functions in this convex combination is an element of 
the set C^^^^y Hence, by (|4.19p . this completes the proof of Theorem [T51 

As in |BGMSG09| we can establish from the above observations strong conver- 
gence of the symmetric velocity gradient and the stress on the subsets 

no ■■= {xen-. \f{(Ti,Si) e A{x) with {cri - S* {x, Du{x))) : {Si-Du{x)) =0 

implies that ^1 = Du{x)}, and 

Qs ■■= {xen-. \f{(Ti,Si) e A{x) with {cri - S* (x, Du{x))) : {Si-Du{x)) = 

implies that cri = S* {x, Du{x))} , 

respectively. Since the proof is identical to the proof of [BGMSG09I Lemma 5.2] we 
omit it here and we only state the result. 

Corollary 22. Assume the conditions of Theorem \18\ and let r' e (l,oo) be such 
that ^ + p- = 1. Then, for all 1 ^ s < r and 1 ^ s' < r' , we have that, as n ^ ao, 

DU"^^u strongly m L'inDf""^, 

^ S*{-,Du) strongly m L^'insY""^- 
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4.3 . Proof of Lemma [201 The proof of this assertion is motivated by the proof 
of |BGMSG09| Lemma 4.6]. However, since we are approximating problem (|1.2p 
with finite element functions here, we need to use the discrete Lipschitz truncation 
from gS31 

Let us define the auxiliary function 

Unix) := DU"ix)) - S\x, Du{x))) : (D(7"(a;) - Du{x)) (4.20) 

and observe that 



Ifln — bn\ dx ^ as n ^ co. (4-21) 
Indeed, thanks to (|2.5p and the properties of r]^, we have that 

f\a„-W,\dx=f f (S'^ix^O-S^ix^Du)) ■.{DU"-Du)diy^{C) 



da; 
dx 



- [ {S*{x,C)-S*{x,Du)) ■.{C-Du)d,y:{C) 
= / / (5*(a;,C)-5*(x,£)«)) :(r>t/"-C)di.:(C) 

^ f f |5*(x,0-5*(a-,r>tt)||r>t/"-C|di/^(C)dx 

in Jm.^^^ 

s; - /" sup \S*{xX)- S*{x,DU''')\dx -. 

Consequently, in order to prove that 6„ — > in measure it suffices to prove that 
a„ ^ in measure. We shall establish the second claim in several steps. 

Step 1: First we introduce some preliminary facts concerning discrete Lipschitz 
truncations. For convenience we use the notation 

E" := n^i^{U" -u) = U" - U^i^u e V" 

and let {E"'-'}nj^fi cz V" be the sequence of Lipschitz-truncated finite element 
functions according to Corollary [T71 Recall from Lemma [TOl that E"^ weakly in 
W^'^'iHy, i. e., we are exactly in the situation of Corollary[T71 Although E"- e Y^w 
i.e., it is discretely divergence-free, this does not necessarily imply that E"''-' e V^iv 
and thus we need to modify E"'^ in order to be able to use it as a test function in 
(|3.17p . Recalling Corollary [5] we define 

^n,j ._ 55n(^j^ _gnj-) ^ y„^ (4.22a) 

The 'corrected' function 

„ ^nj g yn^ (4.22b) 

is then discretely divergence-free. We need to control the correction in a norm. To 
this end we recall from Section [XTI that Q" = spanjQ", . . . , Qtv,,} ^ certain lo- 
cally supported basis. Then, thanks to properties of the discrete Bogovskii operator 
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and Corollary [9l we have that 
/3^||*"'J||i,^ s; sup ^ „^„ = sup 

I 

sup ( 2 



<div £;"^^ Q)^ (div.E;"'J -div£;", g)^ 

Q\\r' QeQ" IIQIIr' 

<div£;"'^ -div£;", p.QfX 



Q=Ili-l PiQi \ suppQf 



<div£;"'^ -div£;", p.grX/ 



/ <div£;"•^p,QfX 



■ III]suppQ7n{£;"'^#£;"}/0Pi'3" 

^ II div£; •■'xa" ||r sup ' 



^c||div£;"^^Xo" . ||.=Sc||V£;"'^xo" . I 
where Yon is the characteristic function of the set 



"£„.,^£„} := U jr^B I £^ e such that E cz {E'''^ ^ E' 

Note that in the penultimate step of the above estimate we have used norm equiv- 
alence on the reference space Pq from (j3.1bp . In particular, we see by means of 
standard scaling arguments that for Q = PiQ? the norms 

(SiP.r'llQrii;;) and Q-ligil.' 

are equivalent with constants depending on the shape- regularity of Qn and Pq. This 
directly implies the desired estimate. 

Observe that |ri^| c |£^| for all E e ^/„, n e N, with a shape-dependent constant 



c > 0; hence, 
Corollary [Tfe) I that 



$1" 



< c|{£;""'' 7^ and it follows from Theorem [12] and 



/3,||*"'^||i,,^c||A„,,xn" „, ||.<c2-J"/'-||V£;"||,. (4.23) 

Moreover, we have from Corollary [T7] and the continuity properties of *8" (see 
Corollary [9]) that 

^"J^ i^n.j _ Q weakly in Wl'^^f for aU s e [1, cx)), (4.24a) 

ipn,3 _^ Q strongly in L'{nf for aU s 6 [1, (»), (4.24b) 



as n ^ CX). 

Step 2: We claim that 



limsup / |a„| dx ^ c2 ^Z'', 



with a constant c > independent of j. To see this we first observe that |a„| = 
a„ + 2a~ with the usual notation a~{x) = max{— a„(x), 0}, x e Q. Therefore, we 
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have that 



hm sup / I a„ I da; ^ hm sup / a„ da; 

(4.25) 



+ 2 hmsup / a„ da;. 

n— >cc J {E'^=E^-i] 

We bound the two terms on the right-hand side separately. As a consequence of 
and the fact that 6„(a;) ^ for a.e. a; e 17 (cf. (|4.12p ) it follows that 



/ a„ d.T ^ / a„ d.T ^ / |a„ — bn\ dx 0, 

J{E"=E"-^ Jn Jn 



as n ^ CO. 



The bound on the first term on the right-hand side of (|4.25p is more involved. In 
particular, recalling the definitions (|4.20p and (|4.22p we have that 



a„ dx 



= / DJ7"(a;)) - S*{x, Du{x))) : {DW^^^^uix) - Du{x)) dx 

J{E" = E"-^ 

+ [ S^i-^DV) : D*"^^ dx-+ [ S^i-.DV) : D*"-^ dx 
Jn Jn 

- f S*{-,Du) : DE"^'dx 
Jn 

+ f {S*i;Du) - S''{;DU'')) : DE''^' dx 

J{E"7^E"'j} 

= r + iv'-^ + iir-^ + IV" + v"^^ 

Thanks to p.4p and (|4.2p we have that, as n oo, 

\r\^ f \S"{x,DU"{x)) - S*{x,Du{x))\\DU2iM^) - Du{x)\ dx 

J {E'^^E"'^] 

< \\ff\-,DU'\-)) - S''{-,Du{-)\\r.\\DlVl^,u-Du\\r ^ 0. 

In order to estimate II"'-' we recall that e V^\^ is discretely divergence-free, 
and we can therefore use it as a test function in p.l7p to deduce that 

ir^ = -B[U'\ [/", *"'^'] + </, *"'^'>j2 ^0 as n ^ cx). 

Indeed, the second term vanishes thanks to (j4.24ap . The first term vanishes by 
arguing as in (|4.7p — observe that for (|4.8p the weak convergence (j4.24ap of 
is suflacient. The term III"'-' can be bounded by means of (j4.23p : in particular, 

lim sup |lir'^'| lim sup||S'(-,£)(7")||r'||D*"'^||r ^ c2-^>, 

where we have used (|4.2p . Corollary [T7I implies that 

lim IV"'^ = 0. 

Finally, by (|4.2p and Corollary [T71 we have that 
lim sup |V"'^| ^ lim sup(||5*(-,Dm)||,- + \\S'\-, DU'')\\r')\\DE'''^ XiE-^E-.^Wr 

2-^/''. 

In view of (|4.25p this completes Step 2. 
Step 3: We prove, for any e (0, 1), that 

lim / \an\^ dx = 0. 



n— >G0 
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Using Holder's inequality, we easily obtain that 

|a„|'' dx = / |a„|''dx+ / |a„|'' da; 




Thanks to (|4.2p . we have that {J^ |a„| dx)^ is bounded uniformly in n and by 
Corollary [17] we have that 

II 



\r 923 r ' 

n.j 



where we have used that {£J"}„gN is bounded in Wq'^{Q,)'^ according to (|4.2p and 
Assumption [HI Consequently, from Step 2 we deduce that 



lim sup f \anf dxfic ^-^^l"- + 



22'-r(l-i;) ■ 



The left-hand side is independent of j and we can thus pass to the limit j — > oo. 
This proves the assertion and actually implies that a„ — > in measure as n ^ oo. 
According to (j4.2ip we have that 6„ ^ in measure and thus we have completed 
the proof of Lemma [201 

4.4. Proof of Lemma 1211 The proof of this Lemma is given in |BGMSG09| . In 

order to keep the paper self-contained, we shall reproduce it here. 

The assertion is an immediate consequence of the result on Young measures 
from |Gwi05| stated in Theorem [31 It therefore suffices to check the assumptions 
therein. The first assumption has already been verified in (|4.16p . 

Step 1: We prove that the sequence {j^"}„eN satisfies the tightness condition. 
From the definition of i/" (cf. (|4.14p ') it follows that 

7"(x):= max |C| = inax |G,(C)K inax (|C| + C)|) ■ 

^GSUpp I/" QGSUpp ^i" ^GSUpp /l" 

We deduce from the inclusion supp /i" cz i?x/r,(-DV"(x)) that ||7"||s c for some 
constant c > and s = max{r, r'} > 1. Since 51 is bounded, the sequence is 
uniformly bounded in L^(ri), and for A/ > we have 



|{x-er!: suppz.;^i3M(0)}| = |{xefl:7"(x)>A/}|< l Ax ^ 



M 



This yields the tightness of {z^"}neN and it follows from part (i) of Theorem [31 that 
is a probability measure, i.e., ||i'2^||5r)t(Rtixd^ = 1 for a.e. a; e 51. 
Step 2: We turn to proving (|4.17p . Recalling (|4.15cp . the assertion follows if 

there exists a nonincreasing sequence of measurable subsets {Ei}i^^ with \Ei\ — > 

a.s i CO, such that the function 

h{x, C) ■■= {s{x, C) - S*{x, Du{x))) : {d{x, C) - Du{x)) 

satisfies (|2.6p on A = H\Ei for each i e N. This can be seen as follows. From (j4.2p 
it follows that ||£>L'""||r. -I- ||-DM||r is bounded uniformly in n e N. Consequently, the 
sequence {c„}„eN, defined by 



c„(x) := c( |D[/"(x)r-^ + \Du{x)r' + ^ + i^)) 



n 

X 



|r>J7"(x)| + + iy (4.26) 

n/ 



is bounded in L^(fl), where k e (fl) from (|2.4p and c > is a constant to be 
chosen later. Hence, Chacon's biting lemma (Lemma [5]) implies that there exists a 
nonincreasing sequence of measurable subsets {iJ jieN with \Ei \ ^ as i ^ oo such 
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that {c„}„eN is weakly precompact in L^{il\Ei) for each i e N. Fix i e N and set 
u> := ^l\Ei. Thanks to the de la Vallee-Poussin theorem (see, |Mey66| ), there exists 
a nonnegative increasing convex function (jj : R4. — > ]R-|_ such that 



(t>(t) 

> <» as t ^ CO 

t 



and sup / (/)(|c„|) dx < 00. (4.27) 

rieN J uj 



Let us also define 

H{x, C) := {S*{x, C) - 5*(^, Du{x))) : (C - r>«(x)). 
By a simple substitution in the spirit of (|4.15p it follows that 



sup / / h{xX)diy^{C)dx 



= sup / / H{x,^)d^i:{^)dx 

< sup sup / / <f>{H{x,$))d^i^{Odx. 

t^R<P[t) neN J ui J {^GRi^^ : \H{x,^)\>R] 

V : ^ ' 

Thanks to the properties (j4.27p of 4> the assertion follows once Jji has been shown 
to remain bounded. To this end, we observe that 

Ji^ sup / sup 0(iJ(x, 4)) dx ^ sup / (/)(|c„|)da;, 

n^N J uj ^<=Bx{Du{x)) n^N J uj 

where we have used that we can choose the constant in (|4.26p so that 
H{x,0 < c (lll'-i + \Du{x)r' + k{x)) (III + \Du{x)\) . 



The assertion then follows from (j4.27p . 

Finally the identities for Du and S follow similarly from the representations 
(j4.15ap and (j4.15bp and the uniqueness of the weak limits (j4.3p and (j4.4p . 

Thus we have completed the proof of Lemma 



5. Conclusions 

We have established the convergence of finite element approximations of implic- 
itly constituted power-law-like models for viscous incompressible fluids. A key new 
technical tool in our analysis was a finite element counterpart of the Acerbi-Fusco 
Lipschitz truncation of Sobolev functions, which was used in combination with a va- 
riety of weak compactness techniques, including Chacon's biting lemma and Young 
measures. An interesting direction for future research is the extension of the results 
obtained herein to unsteady implicitly constituted models of incompressible fluids. 
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Appendix: Auxiliary comments on Lipschitz truncation 

Although similar techniques were used in jBDF12) to prove the properties of the 
Lipschitz truncation, we decided to present a complete proof of Theorem [Ij] for the 
following two reasons: 

• In contrast with the Lipschitz truncation in |BDF12) . the Lipschitz truncation 
in (|3.25p preserves boundary values. This requires changes to the proof that are 
not always obvious. 

• The concept of Lipschitz truncation seems to be new to the numerical analysis 
community. For this reason we have aimed to keep the presentation as self- 
contained as possible. 

Recall the notational conventions introduced in Section 13.51 prior to Theorem 1131 
and the definition p.25p of the Lipschitz truncation. We start with some basic 
estimates. 

Lemma 23. Let A > and v e WQ'^{fiy^ and let {vj}j(^jq cz R'^ be defined as in 
(|3.25bp . We then have, for all j e N, that 

(a) /g** dx < c fg** \Vv\ dx *S cMi\/v){y) for all y e Q** ; 

(b) /g** |V7;| dx cA; 

(c) for fc 6 N with Q* n ¥^ 0, we have 



\vj-Vk\^cf \v - Vj\ dx + cf \v-Vk\dx; 
Jo** Jq** 

(d) for keN with Q* n ^ we have \vj -Vk\iic£j A. 
Proof. We extend v by zero outside il. 



(a) This statement is a consequence of Poincare's inequality and the Friedrichs 



inequality. Indeed, recalling (j3.25bp . for Q* cz Q, we have by Poincare's inequality 
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that 



7 ~"T~^ dxikc-f \\Jv\dx ike -f \\Jv\ dx cM{Vv){y) 

for all y e Q**; the constant c depends only on d. 

In the case Q* cj: $1, it follows from the fact that is a Lipschitz domain and 

Q** = ^J^Q*, that |Q**\ri| 5= c|Q*|, with a constant c > depending on VL. 
Hence v is zero on a portion of Q** whose measure is bounded below by a positive 
constant, which depends on the Lipschitz constant of d^. Consequently, we can 
apply Friedrichs' inequality (cf. |MZ97[ Lemma 1.65]) to deduce that 

\Iv\ dx cM{yv){y) Vy e Q**. 

[(b)] It follows from ](W2)] that (61^2^) n {M.\Ux{v)) = {OdQj) n {M(yv) A} 
0; compare with For z e {0dQj) n {M(yv) A} let := OciVd^ej = 

6'ddiam(Q**); then, 6*^2** c: Br.{z). Consequently, 

f \\7v\ dx sS: c f \Vv\ dx sS: c f \\7v\ dx ^ c M{Vv){z) sS: c X. 
JQ** Je^** JbrAz) 





V 


dx c -f 


/ 




' Q** 




jQ 



(c) Observe that Q* n is equivalent to Qj n Qk and hence we 



obtain from ](W3)] and Qf* = ^^Qf^ie N, that 

\Q** n 2**1 (4V^)"'' max{ I 1,12*1} • 
Therefore, there exists a constant c > 0, depending on d, such that 

- t>fc| sg -f It) - Dj"! dx + -/- |D-t>fc|da; 

^Q**nQ** jQ**nQ** 

^c-f \vj — v\ dx + c -f \v — Vk\dx. 
JQ** JqI;* 

\{d)\ The claim is a combination of (c) (a) (b) and ](W3)] □ 

The next result proves that the Lipschitz truncation is a proper Sobolev function. 

Lemma 24. Let X > 0, v e W^'^ifiy'' and let Vx be defined as in (|3.25p . Then, 
vx-v = V^.(«. - v) 6 W^'\Ux{v) n ny. 

Proof. It follows from (|3.25p and properties of the partition of unity {ipj}j^fq that 
Vx — V = J^j^f^tpjivj — v) pointwise on and vx — v = in the complement 

of Ux{v). Moreover, we have that ipj{vj — v) e Wq'^{Ux{v) n fl^. Indeed, for 
2* c n this follows from the fact that ipeC^{Q*). If on the other hand Q* <^ Q, 

then this follows from Vj = and v 6 W^'^ (Ux{v) n fl). We need to show that the 
sum converges in W^ '^ {Ux{v) n flY. Since fl is bounded, it suffices to prove that 
the sum of the gradients converges absolutely in L^{n)'^^'^. We have, pointwise, the 
equality 

2 ^{Mv, -v))=Y, -v)+ ^,(V^, - Vv)), 

where we have used that both sums are just finite sums, since the family 2* is 



locally finite. Every summand in the last sum belongs to L^{Ux{v) n QY^'^. For a 
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finite subset / c N, we liave, thanks to Lemma [23] and ttie locally finite overlaps of 
the Q*, that 

^i -= [ E -v)+ i^jiVvj - Vv)\ dx 




Note that Ujw\/ 2* ^ and A \Ux{v)\ = A|{A/(Vt;) > A}| < c||Vt;|Ui(n) by 

the weak type estimate p.2ip and v e Wq'^{VlY. Thus, Xu grA/S*'^ ^ XUx(v)^ £ 
i^(R''). Therefore, it follows by Xy q* ^ and the Lebesgue's dominated 

convergence theorem that E/ ^ as / ^ N. Hence the sum Yjjm '^{^ji'^i ^ '^)) 
converges absolutely in L^iVlY^'^ , and the claim follows. □ 



Proof of Theorem [7Si We shall consider parts (a) (d) in the statement of the the- 
orem separately. 

The claim directly follows from vx — v & Wq'^{U\ n VlY (see Lemma [ 
We begin by noting that 



(a) 



(b) 



By Jensen's inequality and the local finiteness of the Q** we then deduce that 
Xu.(.) \vx\' dx < 2 |Q*| ( / |.| dx) ^ 



dx ^ c I \v\ dx, 

n 



for s e [l,oo), which then proves [(b)] for s e [l, oo) using also that vx = v outside 
of Ux(v). The case s = co follows by obvious modifications of the argument. 



(c)]We define := {/s e N : Q* n Q* 0}. Then, on every Q* we have that 
Vt^A = v( 2 ^kVk) = v( 2 ^k{vk - Vj)) = 2 (Vi/;fc)(^;fc - -Uj), 

where we have used that Xifce/ V'fe = 1 on Q*. By Lemma [53] we thus obtain 

Xu>.i.) \^vx\ < c ;X Xq* 2 / ^^^x^ dx < c ;X Xq* 2 / IV^I dx. 

The inequality in part (c) now follows by arguing as in part [(b)| 

p)] It follows fr om the final chain of inequalities in the proof of part [(c)[ above. 
Lemma [23] (b) and the local finiteness of the Q** that 

Xu^iv) |V^;a| cA. 

Since vx = v on 7ix{v), we get the first part of the claim 

\Vvx \ ^ cXxux{v)nn + \Vv\ xu^iv)- 

Recah that Hx{v) = {W^\n) u {M(Vt>) A}. Now, = on W^\VL and |Vt!| 
A/(Vt>) proves that |Vt)| X'Ha(i') ^ ^- This proves the second part of the claim. □ 

The following theorem is the analogue of Corollary [T7] for Sobolev functions. 
Similar results can be found in |DMS08) and |BDF12] . 
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Corollary 25. Let 1 < s < ao and let {e"}„Gti c Wg'''{n)'^ be a sequence, which 
converges to zero weakly in W^'^ (flY , as n — > cX). 

Then, there exists a sequence {A„j}„jg N c M with 2^' <: A„,j ^ 2^'^ "i such 
that the Lipschitz truncations e"'^ := e"^ ^, n,j e N, have the following properties: 

(a) e"'J 6 Wo'"-^\ny and e"'^ = e" on 'Ha„,,; 

(b) ||Ve"'J||oo*ScA„,,; 

(c) e"'J ^ in L=^(ri)'^ as n ^ <»; 

(d) Ve"J — * m L=°(r2)'*^'^ as 7i ^ cx); 

(e) For alln,j e N we have \\^n,jXUx(e")\\s ^ c2~s || Ve"||s, with a constant c > 
depending on s. 

Proof. The assertions follow by adopting the proof of Corollary [T71 □ 
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